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In this contribution 1 review lattice formula- 
tions based on so-called twisted mass (tm) Wil- 
son quarks, focusing on the general properties and 
problems that in my opinion are most relevant for 
on-going studies and applications to physics ^ . 

The simplest of these formulations, which cor- 
responds to QCD with one mass degenerate quark 
doublet, was introduced in ref. as a way to get 
rid of the spurious quark zero modes that plague 
standard Wilson fermions and is referred to as 
twisted mass lattice QCD (tmLQCD). 

1. Nf = 2 twisted mass lattice QCD 

Following refs. |1I2I3| . a sensible lattice formu- 
lation of QCD with an SUf{2) flavour doublet 
of mass degenerate quarks is given by the action 
S = Sg[U] + S^p '> [ij, tp, U], where Sg[U] stands for 
any discretization of the YM action and 

+e-^"^="WcrW+m,]v(x), (1) 
7-V^i^7M(V*+V^), (2) 
W„{r) = -a§ K^t^ + ^'^cr(r) . (3) 



*Based on a plenary talk given at Lattice 2004, the XX 
International Symposium on Lattice Field Theory, held on 
June, 21-26, 2004 at Fermilab, Batavia, IL, U.S.A. 
■^The reader is urged to consult the cited papers for more 
details. For several aspects this write-up should be re- 
garded as an introduction (or a comment) to those works. 



Here Afcr denotes the critical mass, niq is the (yet 
unrenormalized) offset quark mass (such that, 
once Mcr has been set to the appropriate value, 
~ niq as ruq 0), while r G [—1,1] and 
CO G [— TT, tt) characterize the specific Wilson- 
type UV regularization. As discussed below (see 
sect. ll.3|l . Mcr is, up to 0(a) uncertainties, inde- 
pendent of Lu. With the whole 0(a) and numeri- 
cal uncertainty on the critical point brought onto 
Mcrir), one can treat lo and (besides r) as 
exactly known parameters. 

For u) = the familiar action for LQCD with two 
mass degenerate Wilson quarks is recovered. As long 
as 7^ and niq ^ 0, no zero modes of the (two- 
flavour) Dirac operator in eq. which I denote by 
Dp, can occur on any gauge configuration, because 21 

Det[I>F] = det[(5^ + ml sin^ lj] , (4) 
(5 = 75[7-V-f W"cr(r)-|-m, cost^] = Q^ (5) 

The property |(1J solves the problems related to excep- 
tional configurations in the quenched approximation 
as well as to MD instabilities due to exceptionally 
small eigenvalues of in HMC-like algorithms for 
unquenched simulations. 

The extension of unquenched simulation algo- 
rithms of the HMC and multiboson type to the ac- 
tion Q -or eq. 11211 below- is straightforward |^. 
Preliminary results from on-going simulations of tm- 
LQCD seem to hint at a numerical cost comparable 
to that of simulations with staggered fermions 0. 

1.1. Symmetries and reflection positivity 

The symmetries of Nf = 2 tmLQCD are dis- 
cussed in detail in refs. |1I3| . Here I recall only 
the main points arising from those analyses. 

The chirally twisted Wilson term in the action Q 
breaks flavour chiral symmetry in such a way that. 
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if ruq = 0, three symmetry generators are preserved 
(see 0, eq. (4.7)) for any ti;-value. Correspondingly, 
there exist three exactly conserved isotriplet lattice 
currents (see eqs. (4.8)-(4.9)). One of them is the 
neutral vector current, while the other two are mix- 
tures of the charged axial and vector currents (they 
get purely axial at w = ±7r/2 and vector at = 0, tt). 

Charge conjugation symmetry remains exact, but, 
owing to the chiral twisting of the Wilson term, all 
single-axis inversions leave the action invariant only 
if combined with the discrete isospin rotations 7i,2 

71,2 : ^{x) iTia'^{x) , ■4){x) -itp{x)Ti,2 ■ (6) 

or, alternatively, a sign change of the twist angle ui. 
In particular this remark holds for the physical parity 
transformation (here xp = (— x, t)) 



Uoix) -> Uo{xp) 
^.^ Uk{x)-*Ul{xp-ak), fc=l,2,3 
ip{x) — > 'yoi>{xp) 

tp{x) tp{xp)"fo . 



(7) 



For later use I note the non-anomalous symmetry |2| 
7^f = 7^5 X , (8) 
where 



lis 



ip{x) 
ip{x) 



■ isi^ix) 



U^,{x) Uli-x - afi) 



tp{x) 

i^ix) 



g3i7r/2 
3i7r/2 



H-x) 
ip{-x) . 



(9) 



(10) 



Correlation functions evaluated with the tm- 
LQCD action enjoy link reflection positivity for 
all values of r e [—1; 1] and ut G [— 7r,7r) 7 , as 
well as site reflection positivity provided |r| = 1 
and |8r -|- 2aMcr{r) + 2am, cos w| > 6 JTO. 

1.2. Renormalizability 

The tniLQCD fcrmionic action is written 
in what is usually called the "physical basis" |3] , 
with niq real (and positive). By a change of basis, 

the action takes the form considered in ref. PP, 



7- V-l- 



V* +mo+ im2-f5T3 xi^) , (12) 



toq = + cos w , 7712 = TOq sin w . (13) 

All the symmetry properties described above or in 
refs. fll3| can be expressed in the quark basis llllt . 
In particular I note the spurionic symmetries of the 
action II12I I 

TZs X {r ^ — r) x (mo —mo) x (m2 -m2) , (14) 



P X (m2 



-m2l 



with TZs form invariant under IIH and 



r 



Uo{x)^Uoixp) 
U,{x)^Ulixp- 
Xix) ^ ^oxixp) 
X{x) -> xixpho 



ak), fc = l,2,3 



(15) 



(16) 



Standard symmetry and power counting arguments 
imply ,1 that (up to terms irrelevant as a ^ 0) the 
most general action for two-flavour tmLQCD is just 
of the form S,[U] + 5^">[x, X, U], with S,[U] cc Xjgl 
and S^p^ given in eq. 11211 . Invariance under V x 7i,2 
(see eqs. l|7| and ||SJ) rules out parity odd pure gauge 
terms (oc tr\FF\ as a 0). While and m2 
need only multiplicative renormalization, mo under- 
goes additive and multiplicative renormalization 0. 

The continuum flavour chiral WTI's of QCD 
with two mass degenerate quarks can be imple- 
mented (up to 0(a)) for all a;- values |ll3j . with 
the renormalized current quark mass given by 



rhq = Zp'^ZM{uj)mq , 
ZmIi-^) = [ZpZgo cos^ uj + sin^ lu] 



1/2 



(17) 



Zp and Zs't arc the (mass-independent scheme) 
renormalization constants of X75''"aX and XX- 

1.3. Critical mass 

For given values of and r, the appropri- 
ate value of Mcr can be determined by adjusting 
mo so as to enforce one of the properties (chi- 
ral WTI's, pions with minimal mass) dictated by 
flavour chiral symmetry ^ . In the quark basis (|ll|l 



■^In perturbation theory all the conditions of this type 
defining Ma lead to a unique result: in fact 0{ap) terms 
can be singled out and removed, since all momentum scales 
p are controllable at fixed Qq. Non-perturbatively this 
is no longer the case, because there is no way of letting 
ahqcD — » at fixed g^. At > diff'erent definitions of 
Mcr in general difi'or by amounts of order aAqcD ■ 
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and for any m2-value, a sensible condition defin- 
ing is given e.g. by (index 6=1,2 only) 

= 0, (18) 



with O a conveniently chosen multilocal operator 
and X 7^ {y, ...}. Owing to the symmetry 1)15(1 . 
numerical estimates of Met corresponding to dif- 
ferent values of m2 differ by ©(amj) from each 
other. Since, by taking e.g. O = {xibnx){v), the 
parameters ruq and uj enter the condition (|18|l 
only through m2, one sees that, up to irrelevant 
0(a) terms, M^-c is independent of and uo. 

In infinite volume, correlation functions and de- 
rived quantities may depend (e.g. if chiral symmetry 
is spontaneously broken) also on the path along which 
the critical point (mo,m2) = (Mcr,0) is approached, 
and thus on uj. However, since renormalization is a lo- 
cal procedure, any possible tj-dependence of whatever 
determination of Mcr based on infinite volume quan- 
tities is necessarily limited to irrelevant 0{aA.QCD) 
contributions. 

Moreover, the combination of the spurionic invari- 
ances 1141 1 and 1151 1 implies that, if mo = A/cr fulfills a 
certain condition defining the critical mass for a given 
r (and m2), then mo = —Mcr satisfies the same con- 
dition for —r (and the same m2)- The critical mass 
counterterm is thus an odd function of r: 



M„(-r) = -M„.(r), 



(19) 



If there exists an interval [si,S2] of mo-values for 
which the condition defining Mcr(r) is satisfied, the 
invariances 11411 and I15II imply that the mo-values 
in the interval [— S2,— si] are solutions of the same 
condition for — r. In finite volume, analiticity in mo 
excludes the possible existence of such intervals of 
solutions. They may however show up, though with 
a width vanishing as a ^ 0, if Mcr is defined through 
some infinite volume quantity, as for instance when 
it is determined by the vanishing of the charged pion 
mass and the "Aoki phase scenario" |H| is realized 

Once a definition of Mci{r) has been chosen for, 
say, r > 0, one must (and, as shown above, always 
can) take Mcr{—r) = — Mcr(r). Otherwise one 
would unnecessarily spoil -by an artificial choice 
of the untwisted mass counterterm- the spurionic 
symmetry TZc, x V x {r ^ —r) x (mo —itiq) 

'*In this scenario the intervals of mo-values with zero 
charged pion mass are expected to have O(a^) widths. 



enjoyed by the lattice theory prior to renormal- 
ization. For these reasons the criticism to eq. H19|l 
raised in ref. [H] is unjustified. 

2. Mass non- degenerate flavours 

For the case oi uj = n/2 (maximal twist), the 
fermionic action of an SUf{2) pair of mass non- 
degenerate quark is conveniently written .10) as 



-ilbTiWcrir) +mq- T^eq 



7 • V - 

iIj{x) , 



(20) 



where to keep the mass term real and flavour diag- 
onal I have used the matrix to split the masses 
of the members of the doublet. Consequently the 
Wilson term has been twisted with the flavour 
matrix ri. Note that mq and eq are both real. 

It has been shown 10^ that the quark mass 
splitting eq is only multiplicatively renormalized 
and the continuum flavour chiral WTFs can be 
implemented up to 0(a), provided that ^ 



1^ ^ = rhq — iq ~ Zp^m„ — ^e 



S '^9 ' 
Zg^eq 



(21) 
(22) 



are identified as the renormalized (current) 
masses of the quarks in the doublet. 

Remarkably, the fermionic determinant ^, 



det[4;f]>det[Q2^ + m2-e2], 



Qer=75[7-V-KWer(r)] =Qt 



(23) 
(24) 



is real and positive, as long as < ^HI- This 
allows for unquenched Monte Carlo simulations, 
for instance by means of algorithms of the multi- 
boson or PHMC type based on some polynomial 



-1/2 



approximation of [^'pnd-^Fnd 

Alternatively, lattice formulations of QCD with 
Nf mass non-degenerate quarks of the (twisted) 
Wilson type can be obtained by taking the 
fermionic action of the form [11112113] : 



qOS 
Op 



[7-V + e-'^=«/PFcr(r/) +mf]qf{x) } , 



(25) 



^Zs is the renormalization constant of X'''aX: = li2,3. 
''Here -Dpnd denotes the Dirac operator (a 2 X 2 matrix in 
flavour space) corresponding to the fermionic action 1201 . 
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with 9 = ^fdf = '' and Wcr{rf) defined as 
in eq. Q. For all flavours /, the critical mass 
Mcr{rf) can be taken independent of 6*/ and is 
given by the same dimensionless function aMcr{r) 
as in the case of (un)twisted mass-degenerate Wil- 
son quarks |13I14| . The renormalized counterpart 
of the bare quark mass nif reads [11113114] 

rhf = Zm{rf)mf , f = u,d, s,c, . . . . (26) 

Choosing 0/ ^ (at least for u and d quarks) 
allows to get rid of spurious quark zero modes, 
while keeping the action diagonal in flavour space. 
However, even ii 9 = 0, for generic values of 
the TO/'s the lattice fermionic determinant corre- 
sponding to the action (|25|l is complex |11I12I14] . 
Actions of the type can thus be used in gen- 
eral only for valence quarks (see sect-EJ. Partic- 
ular cases with positive quark determinant exist: 
e.g. the Nf = 4 case with m„ = nid, arbitrary nis 
and rrid, 9u = -9d = 7r/2 and 0^ = = [H] , 

3. Continuum and chiral limits 

In tmLQCD the relative magnitude of the lat- 
tice cutoff effects arising from the chiral violating 
action terms may change significantly as a func- 
tion of the quark mass. This remark is partic- 
ularly relevant in the regime where chiral sym- 
metry is spontaneously broken. In this situation, 
in the continuum theory the chiral phase of the 
vacuum is driven by the phase of the quark mass 
term. The same must be true on the lattice, thus 
ideally the continuum limit should be taken -at 
non-zero rhq- before letting rhq 0. 

3.1. 0(a) improvement 

From an analysis a la Symanzik 15 of the lead- 
ing cutoff effects in tmLQCD (eq. (^) and by ex- 
ploiting the invariances |(Hl and ^ 

■Rf = 7^5 X (r ^ ~r) x {ruq -lUq) , (27) 

it follows 3; that the expectation values of multi- 
local, gauge invariant, multiplicatively renormal- 
izable (m.r.) operators O satisfy 

'^The condition S = is necessary for the corresponding 
continuum limit theory be QCD with vanishing "6-term" . 
*The spurionic symmetry T^-g*" holds for generic u) and fol- 
lows directly from the symmetry 1141 using 1131 and 1191 . 



= 2Cg(r;c.)(0(a:i,...,rc„))^-t) +0(a2). (28) 

0(a) effects cancel ^ in the average of results ob- 
tained with opposite values of r (Wilson average). 

In case of maximal twist, ld — ±7r/2, one can 
obtain 0(a) improved results even from one single 
simulation at a given r-value. In fact invariance 
of the action under V x {r ~r) ^'^ implies 

(o(-i>--o)[^^/4)-r/o(o(xf,...))[j:/;;, (29) 

with rjo the formal parity of O, from which the 
second term in the l.h.s. of H28|) can be obtained. 

As eq. (|28ll holds for arbitrary space-time sep- 
arations, all the quantities derived from m.r. lat- 
tice correlators are free from 0(a) cutoff effects, 
once the Wilson average, or - if w = ±tt/2 - the 
average over opposite values of all external three- 
momenta (see eq. ((23): has been taken. 

The extension to mass non-degenerate quarks 
is straightforward, since the symmetries ^ and 
(with obvious modifications) (|27(l remain valid. 

At variance with the familiar Symanzik's pro- 
gramme for 0(a) improvement (151 . this method 
does not require the addition (and the determi- 
nation) of action and operator counterterms. 

The statement that a certain lattice quantity is 
0(a) improved simply means that in the renormal- 
ized theory it approaches its continuum limit with a 
rate (asymptotically) quadratic in a. Moreover, as far 
as one is interested in continuum QCD with massive 
pions, the renormalized quark mass rhq oc mq must 
have a non-zero limit as a ^ Under this con- 
dition 0(a) improvement holds irrespectively of the 
value of rhq . The contrary result of ref . comes from 
the fact that there rhq is allowed to be an 0(a), or 
0{a^), quantity. The question of the magnitude of 
residual cutoff effects, even if parametrically O(a^), 
when rhq is numerically small is however important 
in practice (e.g. for extrapolations of results to a = 0). 

Large cutoff effects can arise when working on 
a coarse lattice at small TOg-values, just because 
both the quark mass term proportional to niq and 
the Wilson term (which has a different chiral ori- 
entation if u; 7^ 0) contribute to the breaking of 

®The same is true also for 0(a^*"'"'^) effects (fc integer). 

This symmetry follows directly from eqs. and IT9t . 
^^If a — > at rhq = and fixed physical volume, no 
spontaneous breaking of chiral symmetry occurs and chiral 
breaking cutoff effects are expected to be harmless 
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chirality (and in particular to the chiral phase of 
the vacuum). In general, to avoid large lattice ar- 
tifacts one should work with parameters such that 
qAq^j^ <C ruq, but for 0(a) improved quantities 
it is conceivable that in several cases the rela- 
tive magnitude of the dominant (as a 0) cutoff 
effects is just 0(a^AQQj3/mq). Scaling tests are 
thus important to assess the magnitude of cutoff 
effects on various observables as a function of rhq. 

3.2. ChPT analyses and phase diagram 

Lattice chiral perturbation theory (ChPT) is an 
expansion in powers of the quark mass (or external 
momenta) and the lattice spacing and provides an ex- 
plicit representation of (lattice estimates of) physical 
observables in the Goldstone boson sector, as well as 
of the chiral phase of the lattice vacuum, in terms 
of a few low energy constants to be determined |16| . 
Such a representation of observables, though usually 
limited to 0(m*), can be very useful in extracting 
physical information from simulation data. When ap- 
plied to tmLQCD [HI, lattice ChPT involves both 
the quark mass rhq and the (rescaled) twist angle 
to — taji~^ [Z go Zp^ tanoj], or equivalently 

nil = rhq cos Lj , m2= rhq sinuj. (30) 

Here I shortly summarize the main results obtained 
to 0{mt), including 0(a) [HI and 0{a'^) 
artifacts, for tmLQCD in different regimes. 

Regime niq ^ 0(a): The chiral phase of the 
vacuum is determined by the term cx rhq, up to 
negligible higher order corrections. 0(a) cutoff 
effects on pion masses and decay constants cancel 
in Wilson averages at generic uj and are automat- 
ically absent at a) = w = ±7r/2. The pion mass 
splitting is O(a^), as expected. Lattice artifacts 
of order a? /rhq however arise from 0(a) terms in 
the chiral Lagrangian, e.g. in the pion masses PT)] . 

Regime rhq ^ 0(a): Once the 0(a) correc- 
tions to the untwisted quark mass (related to the 
intrinsic non-perturbative uncertainty on Mcr) 
are taken into account, the chiral phase of the 
vacuum is still determined by the term oc faq 
(up to small corrections). The concept of 0(a) 
improvement is of course no longer applicable. 

Regime rhq ^ O(a^): After taking into ac- 
count the 0(a) corrections to the untwisted quark 

Dependence on u) disappears as a —> at fixed rhq. 



mass, one finds that O(a^) and 0{rhq) terms in 
the chiral Lagrangian compete with each other in 
determining the chiral phase of the vacuum. De- 
pending on the sign of the coefficient, C2, of the 
term cx o? in the chiral effective potential, one 
finds two possible scenarios which extend to 
twisted mass m2 ^ the "Aoki phase" (c2 > 0) 
and "normal" (c2 < 0) scenarios of ref. [H|. For 
C2 > 0, the Aoki phase transitions are washed out 
into a crossover if rh2 is non-zero. For C2 < 0, the 
first order phase transition in the untwisted mass 
parameter rhi that occurs at (mi, 7712) — (0,0) 
extends itself in the twisted mass plane, end- 
ing with two symmetrical second order points, 
(?7ii,?7i2) = (0, ±?Ti2c), where the neutral pion 
mass vanishes. Note that ?7i2c ~ |c2|a^. The 
value of C2, which depends on (un)quenching, 
and many details of the lattice action, is closely 
related to pion mass splitting: rri^a — oc C2a^, 
with positive proportionality constant. 

4. Numerical investigations 

A first convincing numerical evidence in favour 
of 0(a) improvement via Wilson average and its 
remarkable simplicity in the case of tmLQCD 
with u) = ±7r/2 was obtained in ref. (21]. There 
the scaling behaviour of the vector meson mass 
and the pseudoscalar decay constant was stud- 
ied in quenched tmLQCD with two mass de- 
generate quarks for lattice spacings in the range 
0.068 ^ 0.123 fm. All quantities were expressed 
in units of Sommer's scale rg and, for each value 
of rriQ was set to M^r and m2 was chosen so 
as to obtain romps = 1-79. The results were also 
compared with the corresponding ones for plain 
and clover-improved Wilson quarks. 

An extension of this scaling study to lower quark 
masses (down to vomps ~ 0.6) and few more lat- 
tice spacings is currently in progress |22| and should 
provide useful insights on the size of residual cutoff 
effects as a function of the quark mass. A compara- 
tive study of several hadronic observables for a wide 
range of quark masses, using tm and overlap fermions 
in the quenched approximation ai l3 — 5.85, has also 
been presented |23j . Given the different chiral prop- 

^''Higher order corrections are unlikely to alter the quali- 
tative features of these scenarios I18I19I . 
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erties of tm and overlap fermions, information on the 
continuum limit, coming from scaling tests, such as 
those of ref. |21I22| . may be very beneficial for the 
interpretation of the results. The performance of lin- 
ear solvers for quark propagators in the two different 
lattice formulations has also been studied at /3 — 5.85 
for given volumes and pion masses |24l . 

A nice computation of the pion form factor, F{Q^), 
has been carried out in quenched tmLQCD at to — 
tt/2, 13 — 6.0 and pion masses of about 660 MeV 
and 470 MeV E^. The conserved ( one-point split) 
isotriplet vector current was employed and the 

matrix element 

{n+{pf)\Vo'\n+{p^))=F{Q')[E{pf)+E{p,)], (31) 

with Q = Pf — Pi, was 0(a) improved by averaging 
over opposite values of Pij Results agree with 
meson vector dominance and other 0(a) improved 
computations and cover a wide range of Q^-values. 

The phase diagram of unquenched tmLQCD 
with two mass degenerate quarks (action (|12|) ) 
in the plane (mo, 7712) is currently under study as 
a function of and for different choices of the 
pure gauge action (plaquette and DBW2). Re- 
sults with the plaquette YM action at (3 = 5.2 did 
reveal the presence of metastabilities in several 
observables (plaquette, pion mass, m^^^^) for 
small values of m2 and mo close to its critical 
value (the value where m^*-'^*-^ should vanish) |S] . 

Long living metastable states (associated with val- 
ues ofm^^^*^ of different sign) were identified and the 
metastability in the plaquette was related to that in 
the untwisted condensate (xx)- These findings were 
interpreted as evidence for a segment of first order 
phase transition, i.e. the scenario arising from lattice 
ChPT in the regime rhq ~ O(a^) for C2 < 0. Metasta- 
bilities in unquenched simulations with Wilson-like 
quarks would thus be ultimately related to chiral vio- 
lating cutoff effects and should be absent for m2 7^ 0, 
if C2 > 0, and for \m2\ > m2c ^ |c2|a^, if C2 < 0. 
Whenever C2 < 0, the condition |m2| > m2c can al- 
ways be fulfilled by sufficiently decreasing a, but it 
would of course be desirable to have both C2 small 
in modulus and a good a^-scaling behaviour already 
for a ~ 0.10 -T- 0.15 fm. To what extent this can 
be achieved by simple modifications of the irrelevant 
terms of the lattice action is still an open question. 

-"^^See eq. (12) of ref. [5a] for the definition of m^^^'-^ . 



5. Weak matrix elements and tm quarks 

Although tm Wilson quarks do not preserve full 
chiral symmetry, they can be used to compute 
many weak matrix elements with no or substan- 
tially reduced (with respect to standard Wilson 
fermions) operator mixings. To achieve this re- 
markable result, the key step is to choose the de- 
tails of the UV regularization of the (typically 
four quark) operators on a case by case basis. 
More precisely, one can conveniently employ for 
the valence quarks in the operator of interest 
a flavour diagonal action of the form (|25|l and 
choose the value of the angles Of so as to max- 
imally simplify operator renormalization. If the 
UV regularization chosen for the quarks in the op- 
erator (and the hadron interpolating fields) does 
not admit a positive defined fermionic determi- 
nant, in unquenched studies one can adopt a dif- 
ferent regularization for the sea quarks, e.g. that 
provided by maximally twisted Wilson quarks 
(see eq. (|20() ). In this case, to make contact, as 
a ^ 0, with continuum QCD the masses of the 
sea and valence quarks of the same flavour must 
be matched, for instance by imposing the equality 
of the renormalized current quark masses 

Within the quenched approximation, use of tm 
Wilson quarks and non-perturbative renormaliza- 
tion (in the SF scheme) has recently led to a pre- 
cise computation of Bk, with no mixings and ex- 
trapolation of results to the continuum limit |26| . 

Moreover it has been remarked |I4j that by 
using the clover improved version of the valence 
quark action (P5|l . with Nf = 4, 0„ = —0d = 7r/2 
and either ds = —Oc = tt/2 or ds = Oc = 0, the 
renormalization of K t: matrix elements re- 
quires at most linearly divergent counterterms, 
which can be determined by enforcing parity. 

A rather general strategy to make use of 
maximally twisted Wilson quarks for evaluating 
(un)quenched weak matrix elements with neither 
wrong chirality operator mixings nor 0(a) cutoff 
effects has been presented in ref. Besides us- 
ing different (cleverly chosen) tm Wilson regular- 
izations for sea and valence quarks, the approach 
is based on the remark that renormalizable lat- 
tice models with four sea quarks and a certain 

^^Here I am interested in the theory with massive quarks. 
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number of (possibly replicated) valence quarks, 
plus corresponding ghosts, yield, among others, 
operator matrix elements that coincide -in the 
continuum limit and provided the renormalized 
sea and valence quark masses are appropriately 
matched- with those of Nf — 4 Euclidean QCD. 
The method was illustrated in ref. jl3j by dis- 
cussing the evaluation of Bk as well as K ^ tttt 
and if — i- TT amplitudes. In ref. it has been 
pointed out that it can be extended to static 
quarks and applied e.g. to the computation of Bb- 

6. Conclusions 

Lattice formulations of QCD with chirally 
twisted Wilson quarks provide a framework for 
non-perturbative numerical studies where spuri- 
ous quark zero modes are absent and 0(a) cutoff 
effects, whenever present, can be removed either 
a la Symanzik or by a new simpler method. The 
magnitude of the residual O(a^) scaling violations 
in the small quark mass region is under investi- 
gation. In unqucnchcd simulations a non-trivial 
phase diagram is also expected, and numerical 
studies of it are in progress. Finally, tmLQCD 
can be in various ways extended to mass non- 
degenerate quarks and adapted so as to make pos- 
sible computations of matrix elements of the weak 
effective Hamiltonian with reduced or no mixings. 
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